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WE: 为 提高 空间 堆 中 子 学 模拟 的 可 靠 性 ， 基 于 国际 临界 安全 手册 中 关于 空间 堆 


的 临界 基准 实验 (HMF005 基准 实验 ) ， 采 用 蒙特 卡 罗 程 序 对 HMF005 基准 实验 
中 堆 芯 排 布 进 行 建 模 ， 结 合 不 同 评价 核 数 据 库 ， 计 算 6 种 堆 芯 排 布下 的 kr 计算 
值 与 实验 值 的 偏差 ， 并 对 计算 值 与 实验 值 的 比值 的 平均 值 、 标 准 偏差 、 交 和 均 方 
根 误差 等 进行 深入 分 析 。 研 究 发 现 ，ENDF/B-VIL1、CENDL-3.2、ENDF/B-VIIL0 
和 JEFF-3.3 四 种 评价 核 数 据 库 计算 扒 芯 kg 结果 与 实验 值 偏差 处 于 三 倍 不 确定 度 
之 内 ， 其 中 ENDF/B-VL1 计算 结果 最 接近 实验 值 ，CENDL-3.2 与 实验 偏差 稍 大 。 
CENDL-3.2 库 计算 结果 与 实验 值 偏 差 稍 大 主要 是 由 于 评价 库 中 的 钼 和 皱 核 素 导 
致 ， 钼 的 同位 素 中 "Mo 的 影响 最 大 。 同 时 还 进行 敏感 性 分 析 ，%Mo 和 ?Be 的 总 
上 反应 敏感 性 系数 主要 跟 弹 性 散射 敏感 性 系数 相关 , CENDL-3.2 fv ENDF/B-VIL.1 中 
Be 的 辐射 俘获 敏感 性 差异 较 大 ， 两 者 偏差 可 达 34.9%, 
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随 着 世界 航天 探索 技术 的 发 展 , 目前 , 基于 太阳 能 、 化 学 能 的 空间 电源 技术 能 力 已 经 接 
近 极 限 水 平 , 随 着 空间 活动 不 断 深 入 ,常规 能 源 将 难以 满足 深 空 探测 和 特殊 任务 的 发 展 要 求 ; 
而 核能 因 其 具有 功率 覆盖 范围 广 、 能 量 密度 高 、 受 环境 影响 小 、 寿 命 长 等 特点 ,在 空间 探索 
中 具有 广泛 的 应 用 前 景 。 目 前 , 美国 、 俄 罗斯 和 欧盟 等 国家 都 把 空间 核 动力 技术 视 为 国家 前 
沿 核心 技术 ， 并 积极 围绕 空间 应 用 和 空间 资源 开发 进行 战略 部 署 03。 由 于 利用 数字 反应 堆 
技术 可 以 缩短 反应 堆 研 发 周期 和 降低 研发 成 本 ,目前 各 国都 在 利用 数字 反应 堆 技 术 进 行 空间 
堆 设计 , 评价 核 数据 库 中 数据 的 精度 则 直接 影响 空间 反应 堆 设 计 精 度 , 因此， 开展 评价 核 数 
据 库 应 用 于 空间 核反应 堆 设计 精度 的 评估 对 新 型 空间 堆 的 设计 具有 重要 意义 ,虽然 国 内 外 对 
空间 堆 有 大 量 研究 ， 比 如 前 苏联 在 60 年 代 陆 续 研 制 出 TOPAZ-I 和 TOPAZ-I 等 型 号 的 空间 
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核反应 堆 包 ,美国 在 2018 年 完成 千瓦 级 空间 堆 Kilopower 反应 堆 的 测试 内， 中 国 原 子 能 院 提 
出 SPACE-R ZEHE, 清华 大 学 REAL 团队 开发 的 ACMIR 型 空间 堆 回 , 北京 航天 器 环境 工 
程 研究 所 可 靠 性 与 环境 科学 技术 实验 室 研发 的 ACMIR 型 空间 堆 门 ， 南 京 航空 航天 大 学 的 李 
婷 等 人 提出 的 熔 盐 冷却 空间 堆 概 念 四 等 。 除 了 之 前 的 实验 研究 ， 目 前 很 多 的 概念 设计 大 多 直 
接 使 用 评价 核 数据 库 进 行 空间 堆 反 应 堆 物 理 计算 ,很 少 对 评价 核 数据 库 在 空间 堆 设 计 中 的 可 
靠 性 进行 评价 。 因此 开展 空间 堆 核 数据 宏观 检验 研究 对 空间 反应 堆 的 设计 非常 重要 。 临界 基 
准 实验 通常 被 用 于 评价 物理 计算 模型 、 设 计 方法 以 及 所 用 的 核 数 据 库 的 准确 性 ; 基于 实验 精 
确 测量 反应 堆 的 re 可 实现 对 堆 臣 关键 材料 核 数 据 的 宏观 基准 检验 外。 为 了 对 空间 堆 反 应 堆 
物理 计算 使 用 的 评价 核 数 据 库 进行 临界 验证 ， 本 文通 过 选取 国际 临界 安全 手册 ( ICSBEP ) 
中 HMF005 HEERO, RH MCNP 程序 以 及 不 同 评价 核 数据 库 ， 计 算出 与 实验 对 应 的 堆 世 
kr， 开展 空间 堆 相 关 的 核 数据 入 堆 宏 观 检 验 ， 通 过 分 析 不 同 核 数据 库 ke 计算 值 与 实验 值 差 
值 的 大 小 ， 检 验 目前 评价 核 数据 库 的 可 靠 性 。 利 用 了 MCNP 程序 计算 钼 和 钙 的 敏感 性 系数 
和 不 同 堆 芯 模型 下 的 中 子 通 量 密度 。MCNP 使 用 了 伴随 加 权 的 线性 扰动 理论 计算 敏感 性 系 
1 临 青 基准 实验 

HMF005 系列 实验 为 1987 年 俄罗斯 奥 布 守 斯 克 物 理 和 动力 工程 研究 所 进行 的 临界 实验 。 
实验 按照 堆 芯 材料 排 布 不 同 , 可 以 分 为 六 种 堆 必 。 六 种 反应 堆 堆 世 都 采用 高 浓缩 铀 金属 作为 
燃料 ， 由 于 铀 合金 具有 高 反射 率 、 耐 高 温 、 耐 腐蚀 等 优点 ， 部 分 反应 堆 设 计 采 用 金属 钼 反射 
JZ RTE, 部 分 反应 堆 设 计 也 采用 争 作 为 反射 层 。 以 上 反应 堆 设 计 也 为 检验 钼 和 争 的 评价 
中 子 截面 数据 提供 很 好 的 实验 数据 HMF005 系列 基准 题 包括 六 个 实验 , 本 文 使 用 蒙特 卡 罗 
程序 对 六 个 实验 模型 进行 精细 化 建 模 。 表 1 中 给 出 组 成 实验 中 用 到 各 类 棱 束 的 基本 组 成 材 
料 以 及 每 种 材料 的 原子 密度 , 实验 中 共 涉 及 到 5 种 材料 , 分 别 是 铀 块 、 反 射 层 ( 底部 和 顶部 )、 
钼 棒 、 钙 棒 、 和 空 管 。 这 些 材料 再 组 成 U3 HR. U6 棒 、U8 FE, UII HE, THEE. BOR Ta 
管 等 不 同类 型 的 燃料 或 者 其 它 棒 。 基 准 实验 通 过 排 布 和 组 装 黎 料 棒 、 反 射 层 棒 和 空 管 组 建成 
六 个 堆 芯 实验 模型 ， 具 体 堆 芯 排 布 和 堆 芯 组 成 如 下 图 1- 图 6 所 示 ， 尺 寸 单位 为 厘米 。 

表 1 各 区 域 材料 原子 密度 (1/barn'cm) 


BR 铀 块 顶部 和 低 部 反射 层 钼 棒 pi 空 管 
a 3.3109x104 
SU 3.2968x102 
as 1.4590x 10^ 


228Ly 3.1467x10? 


5.3547x10? 5.3547x102 


3.0875x10? 


Molybdenum 


1.0962x 107 


Beryllium 


5.3293x10? 3.2273x10? 


3.2273x10? 5.3293x10? 


3.2912x10? 


Iron 


8.4596x10^ 


1.3969x10? 


8.4596x10^ 8.4596x10^ 1.3969x10? 


Chromium 


7.5543x 104 4.5748x104 


4.5748x 104 4.5748x 104 7.5543x104 


Nickel 


9.2673x 10% 5.6121x105 


5.6121x105 5.6121x105 9.2673x105 


Manganese 


7.5532x105 4.5741x105 


4.5741x105 4.5741x105 7.5532x10° 


Silicon 


4.2090x 10° 2.5489x 10° 


2.5489 105 2.5489x 10° 4.2090x 10° 


Titanium 


3.7090x105 2.2461x105 


1.0624x10? 2.2461x105 3.7090* 10° 


Carbon 


3.0863x10? 


Aluminum 
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图 1 HMF005_01 堆 芯 排 布 ( 左 ) 和 堆 芯 组 成 CB) 
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HMFO005. 02 堆 世 排 布 ( 左 ) 和 堆 芯 组 成 ( 右 ) 
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图 3 HMF005_03 HEHE (Æ) 和 堆 芯 组 成 〈 右 ) 
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4 HMF005_04 SHEAR (Æ ) 和 推 芯 组 成 ( 右 ) 
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图 5 HMFO005. 05 HERES (Ze) 和 堆 芯 组 成 ( 右 ) 
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图 6 HMFO005. 06 堆 芯 排 布 〈 左 ) 和 扒 芯 组 成 ( f ) 


本 文 首 先 基于 评价 核 数据 库 ENDF/B-VI.1!! , CENDL-32U2 | ENDF/B-VII.0!3! 和 
JEFF-3.304 利用 评价 核 数 据 库 处 理 软件 NJOY2016 制作 了 MCNP 程序 使 用 的 连 


据 库 ， 再 使 用 MCNP 程序 对 有 效 增殖 系数 Cs 以 及 相关 核 素 的 敏感 性 


点 截面 数 


4 


续 


系数 进行 计算 。 


准 偏差 .计算 公式 如 下 : 


AN 


评价 参数 选用 计算 值 与 实验 值 的 比值 的 平均 值 以 及 其 


<C/E>=— - YE (1) 
i=l 
ac = Lt CIE <C/E>)’ (2) 
n-— i-l 
Hop, C WE , 分 别 代 表 对 应 实验 i 的 ks 计算 值 和 实验 值 。 同 时 ，X2 和 均 方 根 误差 


(RMSE ) 也 用 于 评估 计算 结果 和 实验 


RB AE, 表示 计算 值 与 实验 值 的 偏 
模拟 计算 时 设 
代 ， 非 活跃 代为 100 代 。 
3 计算 结果 分 析 
3.1 有 效 增殖 系数 


0. 990 F 


0. 985 


和 人 置 的 温度 为 300K。 在 临界 计算 中 ， 设 置 每 代 10000 个 粒子 ， 


吉 果 之 间 的 一 致 性 ， 具 体 表示 如 下 : 


ICO AE 
= i i 3 
ONE ) (3) 
RMSE = Ey Qn? (4) 
n 
AE, = C, — E, (5) 


差 。 基 准 实验 是 在 300K 温度 下 进行 , 故 本 文 使 用 MCNP 
总 代数 为 1100 


ICSBEP-2014 ”实验 值 
ENDF/B- vi 0 n 算 信 
CENDL-3 KACI 
ENDF/B-VILI inet 
JEFF-3.3 算 值 


e4pen 


图 7 24 


ak 2 HMF005 基准 题 计算 结果 的 x 


F 不 同 评价 核 数据 库 计 算 值 与 实验 值 对 比 


2 , RMSE, <C/E>-1 以 及 A (C/E) 


ENDF/B-Ml.1 CENDL-3.2 ENDF/B- VII.O JEFF-3.3 
x 2 0.686553498 3.945609568 0.972362397 1.850668724 
RMSE 298.290686 715.0881065 354.9903755 489.7414284 
C/E -1 -0.001903728 -0.006670100 -0.002709571 -0.004204959 
A(C/E) 0.002517317 0.005940129 0.002666283 0.002753523 


图 7 给 出 基于 不 同 评价 核 数 据 库 采用 MCNP 程序 对 六 个 实验 计算 的 kr 结果 ， 其 中 基 
准 题 实验 结果 的 不 确定 度 为 360 pem， 六 个 基准 题 蒙 卡 计算 结果 的 统计 误差 最 大 不 超过 21 
pem. PRÉS 7 可 见 ， 部 分 评价 核 数据 库 的 计算 结果 在 1 倍 置 信 区 间 之 内 ,所 有 计算 结果 都 处 
于 3 倍 置信 区 间 之 内 。 根 据 第 二 节 中 的 分 析 方 法 ， 将 采用 不 同 评价 核 数据 库 计 算得 到 的 6 
个 堆 芯 的 kr 结果 与 基准 题 实验 值 进行 处 理 ， 结 果 如 表 2 所 示 。 由 表 2 可 见 ， 采 用 
ENDF/B- VI.1 库 和 ENDF/B-VILO 库 的 计算 结果 与 实验 值 结果 符合 较 好 ,但 基于 CENDL-3.2 
评价 核 数据 库 的 结果 明显 与 其 它 三 个 库 结果 均 有 较 大 差别 ， 均 方 根 误差 (RMSE ) 大 约 是 
ENDF/B-VIILO 计算 结果 的 2 倍 。 由 于 基准 实验 中 大 量 采用 钼 和 争 核 素 做 反射 材料 ， 上 述 差 
异 可 能 是 由 于 反射 层 材料 评价 核 数据 的 不 同 导致。 为 研究 CENDL-3.2 库 结果 偏差 较 大 的 原 
, 将 基准 实验 中 6 种 堆 芯 排 布 的 反射 材料 全 部 去 掉 , 再 使 用 CENDL-3.2 PERI ENDF/B- VI.1 
库 进行 Cr 计算 ， 具 体 实 验 结果 如 表 3 所 示 。 由 表 3 可 见 ， 去 掉 反 射 棒 以 后 ，CENDL-3.2 
库 与 ENDF/B-W.1 库 的 计算 吻合 非常 好 ， 由 此 可 以 推断 ，CENDL-3.2 库 与 其 它 评价 核 数据 
库 的 偏差 主要 自 于 反射 层 材 料 钼 棒 和 钙 棒 评价 核 数据 的 不 同 。 


表 3 HMF005 基准 题 无 反射 层 CENDL-3. 2 和 ENDF/B- VII. 1 k,,, 计算 结果 


01 模型 ”02 模型 03 模型 O4 模型 05 模型 ”06 模型 


CENDL-3.2 0.73984 0.69527 0.66558 0.64135 0.69282 0.71136 
ENDF/B-VI.1 0.73966 ^ 0.69557 0.66535 0.64153 0.69564 0.71098 
Kr 偏差 /pcem — -18 30 -23 18 -18 -38 


在 HMF005 04 HAA 1 根 钼 反射 棒 ， 可 视 HMF005 04 为 全 皱 棒 反 射 。 考 虑 到 
HMF005_01 (全 钼 棒 反 射 ) 和 HMF005 04 (全 争 棒 反射 ) 的 特殊 性 ， 因 此 对 这 两 个 基准 实 
验 进行 深入 研究 。 以 CEDNL-3.2 计算 结果 为 基础 ， 用 ENDF/B-VI.1 库 中 的 钥 和 皱 蔡 换 掉 
CENDL-3.2 库 中 的 铀 和 钙 数 据 ， 然 后 进行 HMF005_01 基准 题 和 HMF005_04 基准 题 计算 ， 
计算 结果 与 实验 值 的 偏差 对 比 结果 如 表 4 所 示 。 由 表 4 可见， 将 HMF005 01 中 的 钥 核 素 
和 HMF005_04 PAPRIKA APIA A ENDF/B- VI.1 库 来 源 后 ， 使 用 CENDL-3.2 库 计算 的 
kf 与 实验 值 的 偏差 明显 偏 小 ， kr 分 别 改 进 643 pem 和 402 pcem。 通过 控制 单一 变量 的 原则 ， 
对 HMF005_01 中 钥 核 素 的 天 然 同 位 素 分 别 进 行 核 数据 库 的 替换 来 研究 钼 各 同位 素 对 ker 
来 的 影响 , 计算 结果 图 8 所 示 。 由 图 8 可 见 , 对 于 HMF005_01 基准 题 , 钼 的 同位 素 中 , ?Mo 
的 影响 最 大 ， 为 301 pom; Mo 的 影响 最 小 ， 其 它 同 位 素 大 约 有 100 pem 的 影响 。 

表 4 计算 值 与 实验 值 ki 的 偏差 (C-E) 对 比分 析 (单位 : 1. 0X10”) 


es CENDL-3.2( 反 射 层 中 的 钼 核 数 据 (HMF005_01) 或 者 锌 数据 
评级 核 数据 库 ENDF/B-VI.1 CENDL-3.2 CHMF005_ 04) 采用 ENDF/B-VII.1 中 的 相应 数据 蔡 换 ) 


HMF005501 457 1038 395 


HME005_04 473 919 517 
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图 8 HMF005_01Mo 各 同位 素 引 起 的 CENDL-3.2 与 ENDF/B-VI.14,, 22 fH 


3. 2 敏感 性 系数 分 析 

针对 HMF005 01 和 HMF005 04 数据 库 替 换 后 带 来 kg 的 差异 ， 对 ”Mo Fil °Be 在 
ENDF/B- VI.1 库 和 CENDL-3.2 库 的 总 反应 、 弹 性 散射 、 总 吸收 和 辐射 俘获 反应 敏感 性 系 
数 进行 研究 。HMF005 01 中 ?Mo 四 种 反应 下 的 总 敏感 性 系数 如 表 5 所 示 ， 对 比 发 现 
CENDL-3.2 和 ENDF/B- WL1 中 ?Mo 的 总 反应 敏感 性 系数 和 总 吸收 截面 的 总 敏感 性 系数 
差异 都 在 1% 以 内 ,弹性 散射 系数 和 辐射 俘获 的 总 敏感 性 系数 偏差 达到 7%。 从 图 9- 图 12 
可 知 ，%Mo 总 反应 敏感 性 系数 主要 由 弹性 散射 敏感 性 系数 决定 ，”2”Mo 核 素 辐射 俘获 和 
吸收 截面 的 敏感 性 系数 较 小 , 在 中 高 能 区 ( 中 子 能 量 大 于 104 ev ) , Mo 在 ENDF/B- VI.1 
和 CENDL-3.2 敏感 性 系数 有 显著 差异 ， 这 可 能 导致 kr 差 值 较 大 的 原因 。 


表 5 CENDL-3. 2 和 ENDF/B- VII. 1"Mo 四 种 反应 道 总 的 敏感 性 系数 


评价 核 数 据 库 ”总 反应 敏感 性 ”弹性 散射 敏感 性 ”吸收 敏感 性 ”辐射 俘获 敏感 性 


CENDL-3.2 3.418 X 102 3.203 X 10? -1.581 X103 . -L506x10? 
ENDF/B-VLI ^ 3421X10? 3.237 X 107 -1.690X103  -1.612X10? 


相对 误差 0.109% 1.04% 6.48% 6.57% 


—— ENDF/B-VII.1 MT=1 一 ENDF/B-VIL1 MT=2 
_— 0010F --- CENDL.32 ,, 0.010 --- CENDL-3.2 
S S 
o o 
RS) 0.008 `S 0.008 
E E 
© 0.006 © 0.006 
oO oO 
20.004 Z 0.004} 
Z B 
= = m 
2 0.002 + 2 0.002 - “1 
o v [ 
Nn 只 ri 
0.000 + 0.000 + == 
4 1 1 1 上 L " 上 L " 1 L 1 1 上 4 上 L 1 
10? 10? 10" 10° 10! 10? 10? 10* 105 10° 107 10? 10? 10 10° 10! 10? 10? 10* 10° 10° 107 
Energy/eV Energy/eV 
92- i E Nia 92. Ry jl mus M 
图 9 "Mo 总 反应 敏感 性 系数 图 10 Mo 弹性 散射 敏感 性 系数 
—— ENDF/B-VII.1 — —— ENDF/B-VII.1 — 
0.00005 MT=-2 0.00005 p MT=102 
--- CENDL-32 --- CENDL-32 
~ 0.00000 — 0.00000 + = 
=i = [=] 
. -0.00005 . 2 -0.00005 } 
9 - g 
iE -0.00010 二 -0.00010 + 
vo 
Š -0.00015 © -0.00015 - 
O O 
>, -0.00020 >, 0.00020 - 
€ B 
-> -0.00025 `= -0.00025 + 
`Ë -0.00030 `Ë -0.00030 + 
= g 
A -0.00035 cà -0.00035 + 
-0.00040 -0.00040 } 
iL £ £ š š £ 3 s t 生生 p p 3. oq 4 p jy 
10? 107 10 10° 40F: 10? 10? 10* 10° 10° 107 10* 10? 10? 10 10° 10! 10? 10? 10* 10° 105 107 
Energy/eV Energy/eV 
92 x = x = 92 =a s +k Fob rb beh. zz oe 
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HMF005_04 中 ”Be 在 四 种 反应 下 的 敏感 性 系数 如 表 6 所 示 ， 对 比 发 现 CENDL-3.2 和 
ENDF/B- VI.1 中 ?Be 的 总 敏感 性 系数 、 弹 性 散射 系数 和 总 吸收 系数 差异 都 在 2% 以 内 ,但 辐 
射 俘获 敏感 性 系数 差别 较 大 ; 图 13- 图 16 分 别 给 出 了 总 反应 截面 、 弹 性 散射 截面 、 吸 收藏 
面 和 辐射 俘获 截面 的 敏感 性 系数 分 布 。 由 图 13- 图 16 ATL, Be 的 总 反应 敏感 性 系数 主要 
跟 弹 性 散射 敏感 性 系数 相关 。CENDL-3.2 和 ENDE/B- VI.I 中 "Be 的 辐射 俘获 敏感 性 差异 较 
大 ， 两 者 偏差 可 达 34.9%。 虽 然 两 个 评价 核 数据 库 "Be 的 辐射 俘获 敏感 性 差异 较 大 , 但 由 于 
辐射 俘获 反应 道 的 敏感 性 系数 较 小 , 占 吸收 截面 的 敏感 性 系数 较 小 , 因此 ,吸收 截面 的 敏感 
性 系数 差别 不 大 。 


表 ó CENDL-3. 2 和 ENDF/B- VII. 1°Be 四 种 反应 的 总 的 敏感 性 系数 


评价 核 数据 库 ”总 反应 敏感 性 弹性 散射 敏感 性 吸收 敏感 性 辐射 俘获 敏感 性 
CENDL-3.2 0.2878 0.267 -4.927 X 10° -2.644 X 104 
ENDF/B-VIl.1 0.293 0.272 4.840 X 103 -1.960 X 104 


相对 误差 1.80% 1.61% 1.81% 34.9% 
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图 15 "Be 总 吸收 敏感 性 系数 16 "Be 辐射 俘获 反应 敏感 性 系数 

3. 3 归 一 化 中 子 通 量 密度 

基于 不 同 评价 核 数 据 库 对 六 个 基准 实验 进行 全 堆 世 的 中 子 通 量 密度 计算 ,计算 结果 如 下 
图 17- 图 22 所 示 。 中 子 通 量 密度 采用 30 群 能 群 结构 ， 具 体 的 能 群 结构 可 以 参考 [10]。 由 图 
17- 图 22 可 见 , 基于 不 同 评价 核 数 据 库 ENDF/B-VI.1 , CENDL-3.2 , ENDF/B-VII.O 和 JEFF-3.3 
计算 的 六 种 堆 芯 下 的 归 一 化 中 子 通 量 密度 趋势 一 致 ， 差 别 较 小 。 由 HMF005_02 基准 实验 到 
HMF005_06 基准 实验 可 以 明显 看 出 反射 层 Be 对 中 子 的 慢 化 作用 ,而 模型 1 相对 于 其 它 基准 
模型 ， 中 子 能 谱 更 硬 。 
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图 17 HMF005_01 归 一 化 中 子 通 量 密度 18 HMF005_02 归 一 化 中 子 通 量 密度 
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图 21 HMF005 05 归 一 化 中 子 通 量 密度 ”图 22 HMF005_06 归 一 化 中 子 通 量 密度 


4 结论 
基于 国际 临界 安全 基准 手册 基准 题 (HMF005 实验 ) , 开展 了 空间 反应 堆 评 价 核 数据 的 

宏观 检验 研究 。 利 用 MCNP 程序 和 四 种 评价 核 数 据 库 ( ENDF/B-W.1 、CENDL-3.2、 
ENDF/B-WL0 和 JEFF-3.3 ) 计算 了 六 种 堆 芯 排 布 空间 反应 堆 基 准 实验 的 kr。 主 要 结论 为 : 

(1 ) 基准 题 的 实验 值 与 蒙 卡 程序 结果 显示 ， 蒙 卡 计算 结果 与 实验 结果 偏差 均 处 于 三 倍 标准 
偏差 之 内 。 其 中 ENDF/B-VIL.1 库 与 实验 值 吻 合 最 好 , CENDL-3.2 库 与 实验 值 相差 偏 大 ; (2) 
CENDL-3.2 库 计 算 结 果 与 实验 值 偏 差 较 大 主要 来 自 于 反射 层 材料 钥 和 皱 核 数 据 。 如 果 采 用 
ENDF/B-VI.1 中 钼 和 外 的 核 数 据 蔡 换 基于 CENDL-3.2 评价 库 的 连续 点 截面 数据 ， 
HMF005 01 ( 全 钥 棒 反射 ) 和 HMF005_04 (全 皱 棒 反 射 ) 的 kr 分 别 可 以 改进 643 pem 和 
402 pcem。 对 钼 的 天 然 同 位 素 单独 研究 后 ， 发 现在 CENDL-3.2 库 和 ENDF/B-W.1 库 中 核 素 
Mo 的 改变 引起 的 kg 变化 最 大 ; (3 ) 通过 对 ENDF/B- VI.1 和 CENDL-3.2 中 ”Mo 和 ?Be 
的 敏感 性 系数 进行 研究 ， 发 现 ”Mo 和 ?Be 的 总 反应 敏感 性 系数 主要 跟 弹性 散射 敏感 性 系数 
相关 ,CENDL-3.2 fil ENDF/B- VII.1 中 *Be 的 辐射 俘获 敏感 性 差异 较 大 ,两 者 偏差 可 达 34.9%。 
从 六 个 基准 实验 的 中 子 能 谱 可 见 ，?Be 对 中 子 的 慢 化 效果 明显 。 
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Study on criticality benchmark testing in spatial reactors 


based on different evaluated nuclear databases 


Chen Xiang! Ma Xubo! Fu Liantaol Ma Xudong! 
(1.North China Electric Power University, School of Nuclear Science and Engineering, Beijing 102206) 


Abstract: In order to improve the reliability of neutronics simulation for space reactors, based on the criticality 
benchmark experiment for space reactors in the International Criticality Safety Handbook (HMF005 benchmark 
experiment), the Monte Carlo program is used to model the core arrangement in the HMF005 benchmark 
experiment, and the deviation of the ke calculated values from the experimental values under six core 
arrangements is calculated by combining the different evaluated nuclear databases and the ratios between the 
calculated and experimental values are analyzed in depth. mean, standard deviation, y? and root mean square error 
are analyzed in depth. It is found that the deviation of the core ke results from the experimental values calculated 
by the four types of evaluation nuclear databases, ENDF/B-VII.1, CENDL-3.2, ENDF/B-VIII.0 and JEFF-3.3, is 
within the triple uncertainty, with the calculated results of ENDF/B-VII.1 being the closest to the experimental 
values, and the deviation of CENDL-3.2 from the experiments being a little bit larger. The deviation of the 
calculated results from the experimental values of the CENDL-3.2 library is also shown in the table below. 
CENDL-3.2 library calculations deviate slightly from the experimental values mainly due to the molybdenum and 
beryllium nuclides in the evaluation library, and the isotope of molybdenum has the largest effect of 2Mo. 
Sensitivity analysis was also carried out, and the total reaction sensitivity coefficients of ?Mo and ?Be were mainly 
related to the elastic scattering sensitivity coefficient, and the difference in radiation capture sensitivity between 
CENDL-3.2 and ?Be in ENDF/B-VII.1 was large, and the deviation between the two could be up to 34.9%. 

Key words: critical calculations; effective multiplication factor; evaluation of the nuclear database; Monte Carlo 


simulation 
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